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The Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) states for two-dimensional s- and d-wave super- 
conductors (s- and d-SCs) are self-consistently studied under an in-plane magnetic held. While the 
stripe solution of the order parameter (OP) is found to have lower free energy in s-SCs, a square 
lattice solution appears to be energetically more favorable in the case of d-SCs. At certain sym- 
metric sites, we find that the features in the local density of states (LDOS) can be ascribed to two 
types of bound states. We also show that the LDOS maps for d-SCs exhibit bias-energy-dependent 
checkerboard patterns. These characteristics can serve as signatures of the FFLO states. 
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The inhomogeneous superconducting state known as 
the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state was 
predicted in the mid-1960's for a superconductor (SC) in 
a strong exchange field, due to ferromagnetically aligned 
impurities . When a high magnetic field (H ) is applied 
parallel to the layers of a quasi-two-dimensional (q-2D) 
SC, the Zeeman effect dominates over the orbital effect. 
Then the FFLO state can also have lower free energy 
than a homogeneous superconducting state 0, • 

A number of q-2D High-T c SCs and q-lD and -2D or- 
ganic SCs have been suggested to have the FFLO phase 
in the low-temperature (T) and high (H) limit 0, H, IS 
01 IB- 111 EH Recent interest is on the q-2D heavy fermion 
compound CeCoIn R , w hich becomes a d-wave SC (d-SC) 
below a T c ~ 2.3K [H El El ■ With a strong H applied 
parallel to its conducting planes, heat capacity measure- 
ments reveal a second order phase transition into a new 
superconducting phase, which strongly suggests that the 
system may have realized the FFLO pha se 0, ■ Vari- 
ous other measurements 0, 0, 0, 0| also support the 
existence of the FFLO state. 

The direct observation of superconducting order pa- 
rameter (OP) modulation is a challenging task. The 
quasiparticle local density of states (LDOS) is a useful 
quantity to measure for this purpose. It could be di- 
rectly detected in scanning tunneling microscopy (STM), 
which measures the local tunneling conductance between 
a normal- metal tip and a SC. Vorontsov et al. |2(1 | 
have calculated the LDOS of 2D d-SC by solving quasi- 
classical Eilenberger equations, but only ID OP modu- 
lations have been considered. References [21] and [22] 
treated the OP as a small parameter near the phase 
boundary, and predicted for d-SC that the energetically 
favored state is a square lattice, but they did not in- 
vestigate the LDOS. In this letter, using a tight-binding 
model, we have self-consistently determined the spatial 
variation of the FFLO-state OP in the low T limit, for 
both 2D s-SCs and d-SCs, with an H applied paral- 
lel to the conducting plane. We have then calculated 



the magnetization density and the LDOS. Our results 
show that for 2D s-SCs (d-SCs), ID-stripes (2D-lattice) 
solutions are more energetically favorable. As a signa- 
ture of ID-stripes solutions for s-SCs, the LDOS spec- 
trum shows two low-energy peaks on top of a low-energy 
bump, which indicate a midgap-states |23]-formed mini- 
band with higher LDOS near the edges. In addition to 
two low-energy peaks arising from similar midgap-states, 
LDOS spectrum for d-SCs shows two additional low- 
energy peaks, which are outside the other two peaks in 
energy, and originate from finite-energy Andreev bound 
states (ABS) localized around saddle points of the OP, 
providing definitive signatures of 2D modulations in the 
OP. Measuring hxed-bias spatial maps at these peak en- 
ergies can confirm this interpretation. 

We begin with a model system on a 2D square lattice 
with a pairing interaction V between two electrons on 
the same site for s-SCs, or on the nearest-neighbor sites 
for d-SCs, which in the mean-field approximation leads 
to the discrete Bogoliubov-de-Gennes (BdG) equations: 
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where the single particle Hamiltonian Hij, a — —Uj — 
(fx + ah)Sij with /i the chemical potential, h the Zee- 
man energy or the exchange interaction, and t the ef- 
fective hopping integral between two nearest-neighbor 
sites. u^g.) are the Bogoliubov quasiparticle ampli- 

tudes on the j-th site. The superconducting OP satis- 
fies the self-consistency condition Ay = ^ J2n( u ?i v j'* + 
u^v? l *)tanh{E n /2k B T). The s-SC OP at site % is A u 
and the d-SC OP at site i is defined as Aj = (A< ^+ es + 
A M _ e:c - A M+ej( - A M _ e J/4, where e XiV denotes the 
unit vector along x, y direction. For a 2D system with 
the magnetic field parallel to the superconducting plane, 
there is no orbital effect. (That is, the Doniach phase 
coupling 13] between neighboring layers of q-2D layered 
superconductors is neglected.) In case there are more 
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FIG. 1: Plots of the OP An as a function of position for an 
s-SC. For (a), V = 2.0, h = 0.25; for (b), V = 2.5, h = 0.4. 
Here the chemical potential \x = —0.4. 



than one solutions for a set of parameters, we compare 
their free energies 25] in order to find the energetically fa- 
vored state. We apply this method in real space to study 
the FFLO state with a real OP. A weak point of our 
approach is that we are unable to obtain solutions with 
periods not commensurate with the lattice size. That 
the lattice size used here can not be too large prevents 
us from obtaining solutions for small h and the complete 
phase diagram. The advantage of our approach is that it 
allows us to obtain the LDOS accurately and the physics 
associated with the spatial structure of the LDOS can be 
understood in terms of quasiparticle ABSs. 

The OP is calculated self-consistently by solving the 
BdG equations iteratively, with the OP initially uniform 
on all sites. The calculation is made on 16 x 16 supercells 
each of which is a 32a x 32a lattice and 8 x 32 supercells 
each of which is a 64a x 16a lattice, with a the lattice con- 
stant. The value of the OP at H = is determined by the 
chemical potential fi and the pairing interaction V. The 
modulation period of the OP is close to JTh ^ 1 [Tij , where 
v j is the Fermi velocity, which is determined by \x. Thus 
one can adjust the period by changing either \i or h. To 
change h, we have to change V as well, because an FFLO 
state exists only for h in an interval [fioH c i, /J,oH C 2] Q, 
where /io is the magnetic moment of an electron, and 
the lower and upper critical fields H c \ and H C 2 are both 
proportional to the superconducting gap at zero H . In 
order to obtain solutions that are commensurate with the 
lattice, we relax the restriction on /i and V from the val- 
ues appropriate to any specific material. Throughout our 
calculation, we let t = 1 and T = 0.001. 

The energetically favored solution in s-SCs always 




FIG. 2: Plots of the OP A, as a functions of position for a 
d-SC. For (a), V = 1.0, h = 0.15 and fj, = -0.4; for (b), 
V = 2.0, h = 0.4 and fj, = -0.4. 



forms ID stripes, which fits closely to the form 
Ao cos(gx), with Ao the maximum value of the OP which 
is not necessarily the value of the OP at zero field. We 
use a 64a x 16a lattice to obtain longer period solutions. 
Figure n shows the OP A« as a function of position for 
two choices of V and h at fi = —0.4 for s-SCs. The ID 
modulation of the OP with a periodicity 64a/3 is plotted 
in Fig. ITfa) with V = 2.0 and h = 0.25. Using V = 2.5, 
we see in Fig. ^b) that the periodicity becomes 16a at 
h = 0.4. We find that the periodicity decreases from 16a 
to 64a/5 as the Zeeman energy increases from h = 0.4 to 
0.45. For V = 2.0 and V = 2.5, the OPs at zero field are 
0.34i and 0.57i respectively. The nodal lines where the 
OP changes sign are along (010) or the y direction. 

The OP of the energetically favored states in d-SCs 
always forms a square lattice, which fits very closely to 
the form of A [cos(gx) + cos(gy)], as shown in Fig. |21 
As a consequence, the nodal lines of the OP are in (110) 
and ( 110) directions. The lattice used for our numerical 
calculations is 32a x 32a. ID modulation solutions for 
d-SCs by using a 64a x 16a lattice are also obtained, 
but they have higher free energy per site than 2D lattice 
solutions on a 32ax 32a lattice. Fig. Eta) shows the profile 
of the OP at V = 1.0. Here, the OP at zero field is 0.057t. 
In order to show the periodic structure of spatial profiles 
more clearly, for this set of parameters, we have plotted 
our results on a 64a x 64a lattice by spatial repetition. 
Fig-HM presents the profile of the OP at V = 2.0. For 
this V, the OP at zero field is 0.23t. 

The magnetization TOj=n;f-n^, with (n^) as the 
number of spin-up (-down) electrons at site i, has mod- 
ulations shown in Fig. [21 The parameters of Fig. [21 (a) 
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FIG. 3: The spatial profile of magnetization m; of the FFLO 
state, (a) is for an s-SC with V = 2.0, h = 0.25 and ^=-0.4. 
(b) is for a d-SC with V = 1.0, h = 0.15 and (i = -0.4. 



are the same as those in Fig. (a). The parameters of 
Fig. |3(b) are the same as those in Fig. 0(a). The period 
of modulation here is one half of that of the OP. The 
magnetization in an s-SC (a d-SC) forms ID stripes (2D 
lattice). The magnitude of the magnetization is highest 
at the nodal lines of the OP. 

The LDOS of spin-up and -down quasi-particles can be 
calculated from 

Pia{E) = Y,\K«\ 2 5{E n — E) + \v&\H(E n + E)] (2) 

n 

where the delta function 5{x) has been approximated by 
e/?r(a; 2 + e 2 ). e is set to 0.01 in our calculation. 

In Fig. EI a), we plot the spin-up LDOS for an s-SC 
at three sites with reference to Fig. ^a). Curve A is 
for a site on the nodal line, curve B is for a site where 
|Ajj| has the maximum absolute value, and curve C is 
for a midpoint between sites in A and B. In the upper- 
right corner of Fig. 0] we have the contour plot of the 
OP for a d-SC with reference to Fig. EI a). Here a is 
the saddle point where two nodal lines intersect, b is the 
mid-site between two neighboring saddle points, c is the 
site where the OP has the maximum magnitude, and d 
and e are respectively the mid-points between a-c and 
b-c. The spin-up LDOS at these symmetric sites are 
shown in Figs.^Ib) and (c). The spin-down LDOS spec- 
tra (not shown) are simply the curves in Fig. 0] shifted 
to the right by 2h. Notice that there is a van Hove peak 
located at E = 0.17 (E = 0.27) in Fig. Ha) [Figs. Efb) 
andQJc)]. The LDOS spectrum is skewed by this van 
Hove peak. These plots reveal that there are two kinds 
of ABSs in the d-wave case: One due to the sign change 
of the OP across the nodal lines, and is similar to the 
midgap states |23{. These states would have essentially 



FIG. 4: The spin-up LDOS spectrum in the FFLO state, (a): 
the LDOS at three different locations (A, B, and C, see text,) 
for an s-SC with V = 2.0, h = 0.25 and = -0.4. In the 
upper-right corner, the contour plot of the OP for a d-SC is 
displayed with five symmetric sites, (b) and (c) : Plots similar 
to (a) for a d-SC with V = 1.0, h = 0.15 and ^ = -0.4 at 
sites a, b, c, d and e (see text). The locations of the LDOS 
peaks are marked by arrows. 



zero energy if not for the Zccman shifts to ^fh for spin-up 
and -down quasiparticles, and broadening into a mini- 
band due to tunneling between the different nodal lines 
in a periodic FFLO state. This feature already exist in 
an s-SC, as shown in Fig. 01a), and is the signature of a 
ID FFLO state. The semiclassical orbit of such a state is 
a round trip along a straight-line segment across a nodal 
line, ending with two Andreev reflections involving op- 
posite signs of the OP on the two sides of the nodal 
line. The same mechanism also gives rise to the inner two 
peaks in Figs.^Jb) and^c) which are for a d-SC. The 
second kind of ABSs is essentially localized at the saddle 
points of the OP. Because the FFLO state in d-SC is a 
2D lattice, the absolute-OP is very much suppressed in 
a cross-shaped region at every saddle point, which yields 
a potential well for a quasiparticle and thus generates 
two finite-energy ABSs (per spin) at every saddle point. 
Zeeman shifts also apply to these states, and communi- 
cation between the saddle points turns these states into 
two narrow mini-bands. The semiclassical orbits of these 
states are round trips along straight line segments across 
the saddle points, ending with two Andreev reflections 
involving only one sign of the OP. Thus these are not 
related to the midgap states, and should not have zero 
energy before Zeeman shift and broadening. The LDOS 
plot at a saddle point [the curve a in Fig. Efb)] exhibits 
two strong peaks corresponding to two ABSs located at 
—h ± eo with eo = 0.10, where —h is the Zeeman shift. 
One can see that curve a is modified by two weak-inner 
peaks originated in the midgap states. Because of non- 
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FIG. 5: The LDOS maps for spin-up quasiparticles at certain 
low-energy-peak biases E for d-SC with V = 1.0, h = 0.15 
and (j, = -0.4. For (a), E = -0.13. For (b), E = -0.25. 

locality of these states, 4 low-energy peaks show up at 
various symmetric points in Figs. 4(b) and 4(c). They 
are located at E = -0.25, -0.17, -0.13 and -0.05. (Re- 
moving the Zeeman shift there would be at E = ±0.10 
and ±0.02, representing two ±E pairs of ABSs of differ- 
ent nature.) 

For an s-SC, the LDOS map for spin-up quasiparticles 
at either of the two low-energy peaks, not presented here, 
forms simple ID stripes and its intensity is highest at the 
nodal lines of the OP, as is expected for zero-energy ABS 
(or midgap states) formed near the nodal lines. 

For a d-SC, the structure of the spin-up LDOS maps 
at the four low-energy-peak biases are richer in physics, 
and are presented in Figs. 0a) and0b). Fig. 0a) shows 
the LDOS map at E = —0.13 (with similar result at 
E = —0.17.) The intensity is highest at points on nodal 
lines halfway between two neighboring saddle points, and 
is lowest at both the saddle points and the absolute- 
OP maxima. The high-intensity spots form a 16a x 16a 
checkerboard pattern with unit vectors along the (100) 
and (010) directions. Figure 0b) is at E = —0.25 (with 
similar result if E = —0.05). The highest intensity is at 
saddle points of the OP, and they form a 16\/2a x 16V2a 
checkerboard pattern but the unit vectors are along (110) 
and (110). Thus the two peaks at E = -0.17 and -0.13 
correspond to zero-energy ABS formed near the nodal 
lines, similar to the case in the s-SC. On the other hand, 
the two peaks at E — —0.25 and —0.05 correspond to 
finite-energy ABSs formed near the saddle points. Note 
that if the STM tip is unpolarized, then the measured 
LDOS should come equally from the spin-up and spin- 
down contributions. Even under this situation our calcu- 
lation indicates that the LDOS maps remain qualitatively 
unchanged as compared to Fig. 

In Summary, using a tight-binding model, we have 
studied, in the low T limit, the FFLO state in a 2D SC 
with a magnetic field applied in plane. The supercon- 
ducting OP is self-consistently determined. The spatial 
profile of the energetically-favored solution is found to 
form ID stripes for an s-SC, and a 2D lattice for a d- 



SC. The spatially- varying magnetization and the quasi- 
particle LDOS spectrum have been calculated. At all 
symmetric sites of the FFLO lattice, we find that the 
low-energy features in the LDOS can be ascribed to two 
types of ABSs, one is related to zero-energy ABSs (or 
midgap states) and the other is not. LDOS maps are 
presented for a d-SC at certain low-energy-peak biases. 
They are shown to have bias-dependent checkerboard 
patterns which are different from the spatial profile of the 
OP. Measuring such maps can clearly reveal the nature 
of the two different types of ABSs. These characteristics 
provide clear signatures of the 2D FFLO state in d-SCs. 
In case the phase of CeCoIn 5 has a competing antiferro- 
magnetic order, we need to use a model suitable for the 
heavy fermion system to reexamine this problem. 
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